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ABSTRACT 


The implementation of a no-drop bomb simulation system 
on a micro-computer is presented. The bomb trajectory algo- 
rithms supplied by Mr. Lee Thomson, China Lake Naval Weapons 
Center, formed a base for the simulation programs. 

The simulation program was written in FORTRAN to test 
the four phases of the system: tracking, position and veloc- 
ity fitting, trajectory simulation, and impact point calcu- 
lation. Translated PL/M programs form a base for the 


micro-computer implementation. 
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I. INTRODUCTION 


The No-Drop Bomb Simulation on a micro-computer, con- 
sists of a series of algorithms which are programmed to sim- 
ulate the fall of a bomb in a practice environment. The 
advantages of such a system are numerous, both in cost and 
flexibility. 

Presently, the Department of the Navy spends in excess 
of $1 million for the purchase of practice bombs for pilot 
training.* These bombs require storage space on the air- 
craft carriers for at sea training. Since these practice 
bombs are live ordnance, there are many restrictions placed 
both on practice ranges and training flight plans. 

A simulation of the bomb drop will eliminate the need 
for these practice bombs, in exchange for a simulation micro- 
computer interfaced with a radar, self-contained in a port- 
able van. The van will then be positioned near a designated 
target location prior to the practice bombings. Since there 
is no live ordnance attached to the aircraft, the targets 
themselves can be more realistic; for example bridges, build- 
ings or industrial sites. The pilot's primary restriction 
in number of runs will be the amount of fuel he carries. He 


may Simulate bombing runs as often as the training session 


1 Based on oral communication with the Naval Test and 
Evaluation Squadron (VX-5), China Lake, Calif. 
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dictates with a variety of bombs. Training session time is 
increased due to the elimination of bomb loadings and mechan- 
ical failures. 

The biggest advantage, however, exists in the actual 
training of the pilot. There are only three pilot controlled 
factors which determine a target hit or miss: dive angle, 
velocity, and altitude. With practice ordnance on board the 
aircraft, an uncontrollable variable is introduced-ballistic 
dispersion. The random ballistic dispersion serves only to 
confuse the pilot. Therefore, the release conditions should 
be the only variables that determine the practice run out- 
come. | 

The simulation process is subdivided into four phases: 
the tracking of the aircraft, determining the aircraft's 
position and release velocity, Simulating the bomb fall, and 
calculating the impact point. (See Figure 1). 

A Nike Hercules radar will continuously track the air- 
craft in a spherical coordinate system. An analog to digital 
converter will provide digital input to the micro-computer. 
The micro-computer will store the current position history 
of the aircraft in three arrays, each containing the sixteen 
most current spherical coordinate values. (Section II-A). 

When the pilot releases the bomb, a tone Signal will be 
sent to the micro-computer. A least squares polynomial fit 
using the current aircraft position history yields the spher- 
ical coordinate release position and velocity. This position 
and velocity is translated into a Cartesian coordinate sys- 


tem. (Section II-B). 
10 
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Using the calculated release conditions, the bomb drop 
trajectory is assumed to be governed by a system of differ- 
ential equations which are solved numerically using the 
Runge-Kutta fourth order integration scheme. The bomb alti- 
tude is continuously compared with the target's altitude for 
impact. (Section II-C). 

The impact coordinates are finally calculated using the 
bomb's down range travel and time ok fail. The velocity of 
the wind is taken into consideration in deriving two sets of 
coordinates for the impact point, one with respect to the 
target's fixed run-in-line and the other with respect to the 
aircraft's run-in-line. (Section II-D). 

The simulation was initially programmed using FORTRAN 
for ease of testing end debugging. Important phases in pro- 
gramming were in calculating the aircraft's release condi- 
tions and simulating the bomb drop trajectory. 

The trajectory phase of the process was tested independ- 
ently for accuracy in results. This made it possible to use 
known release conditions for an output comparison with pub- 
lished results. The routine was subsequently optimized for 
the micro-computer implementation. (Section II-A). The re- 
lease condition curve fitting process was similarly tested 
using actual recorded data as the standard. 

Upon completion of the FORTRAN implementation, the re- 
sultant program was translated into PL/M, a high level lan- 
guage compatible with the selected micro-computer (INTELLEC 


8). Because suitable arithmetic operations were not 
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available, systems implementation involved considerable 
additional effort in the development of a floating point 
package for the arithmetic operations. (Section IV-A). 

Due to time limitations, the project came to an end with 


all the PL/M routines compiled but untested. 
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II. BALLASTIC PROBLEM 


The baiiastics problem is subdivided into four distinct 
phases: tracking, position and velocity fitting, bomb trajec- 
tory Simulation, and deriving the impact point. A radar 
tracking system provides the micro-computer with a digital 
spherical coordinate history of the aircraft's position at a 
predetermined time interval. At the bomb's release point, a 
least squares data fit allows the calculation of the Carte- 
Sian position and velocity components at release. The 
bomb's trajectory towards the target is then calculated yield- 
ing a down range travel and time of fall of the bomb. The 
impact point is then derived with respect to the fixed tar- 
get's run-in-line and the aircraft's run-in-line. The ini- 
tial release conditions and the impact point are then relayed 


to the pilot. 


A. TRACKING PHASE 

The initial phase of the overall simulation involves 
tracking the aircraft until the pilot releases the bomb. The 
tracking phase samples the radar's special spherical coordi- 
nate inputs at a fixed time interval yielding R, 0, and i 
(See Figure 2). 

These input values are stored in three wrap-around arrays 


corresponding to R, 8, and >. In order to minimize the 


Special coordinate system based on the radar's output. 
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Figure 2, Radar Tracking System Coordinates 


arithmetic processing and storage, the values stored in these 
arrays are the actual radar outputs subtracted from the base 
value. Since each of the arrays are cyclical in nature, new 
inputs are read over old input values thus retaining the 
most recent sixteen point position history of the aircraft. 
Rae array consists of sixteen elements and two associated 
base values: Tilda and Bar. Two base values are needed for 
the wrap-around type input. (See Figure 3). 

| When the table is full, a new base value is calculated 
and the old base value is set equal to the current base val- 
ue (R « R). The current base value is set equal to the cur- 
rent radar output value (R < R), and the value stored in the 
table is equal to the current base value minus the current 
radar output value (A Re R). 

Once the pilot releases the bomb, a current spherical co- 

ordinate position history is available with the associated 


base values. 


B. POSITION AND VELOCITY FIT PHASE 

To reduce the number of arithmetic operations, the least 
Squares position and velocity fit was performed on the spher- 
ical coordinates and subsequently translated to Cartesian 
position and velocity components. 

The least squares fit was performed using the second de- 
gree Legendre Discrete Orthogonal Polynomials. The second 
degree polynomial was used to estimate the position and ve- 


locity of the aircraft at the release point. The Legendre 
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Orthogonal polynomials of the second degree were used for a 
sixteen point fit. = 
The orthogonal polynomials have the form: 


9 15 6*) = 1 


G7, 15) Peal Be 
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15 
1-6+* +6 * X + eo 
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a HANS = .GFTGEX + 6EXCK= 1) 
14*15 


The values of these polynomials are based on the number 
of points to be fit. (See Table I for polynomial tabula- 
tions). 

The second degree polynomial which fits the data is 
given by the formula below. For further reference see 
McCalla [1]. 


F(x) = a, + a, (1 = Ue) + a, (1 ZN (greer) g+X_4{x-1)) 
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SECOND DEGREE LEGENDRE DISCRETE POLYNOMIALS EVALUATED FOR A 16 POINT FIT 


Table I. 
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where Track (I) refers to one of the data arrays R, 6, or 9. 
For example, if release occurs at position K in the track 


array,0O is multiplied by Track(K+1). Since the track 


0,0 
array values are differences, the proper base value must be 
applied. Therefore, for release at point K, the a.'s 


(j « 0,1,2) are calculated as follows: 
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It can be seen that the sum of the products using the Track 
. Gifferences minimizes the arithmetic processing. 


To evaluate F(x) at x=15 (release time) we obtain: 
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Smoother position at release. 
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To evaluate F'(x) at x=15 we obtain: 
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The special spherical coordinate positions and velocities 
are transformed to the CarteSian coordinate positions and 
velocities using the following formulas: (See Figure 4). 


X =R * Sin 6 * Cos 6 


Vv 
Yy ent COS iG +) Cosi @ 
Z =R* Sin o 
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thus taking time derivatives: 


xe = Sin @*Cos $*R + R*Cos @*Cos o*6 + R*Sin 9*Sin o*$ 
ve = Cos 6*Cos $*R + R*Sin @*Cos o*6 - R*Cos 9*Sin o*d 
Za = Sin 9*R + R*Cos 0*6 


The reader will note that the Xo yy? and Za position co- 
ondiaates are with respect to the radar's (van) axis. The 
positions are translated by using the van's displacement co- 
ordinates from the target center Oe eae and CANNED) to 


put position coordinates in the target's frame of reference. 
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(See Figure 5). 
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Figure 4, Transformation of Radar Tracking Coordinates 


into Cartesian Coordinates, 


iw) 
to 


syotinvesdD suritoar? ashen te ve £9 Barse Derugd a 


»te®an iPro) malasrrsed ofnt 


- 
. 


Sowely ofeCUTPI00D °C sansTy 


uBsA 
03 300dsez /m uot,]}sod s,ayeaoaty - Ags ee 


e974 
07 yoodsez /m uoy,tsod s,4yeqouty - 7 ' x 4 yy 


UBA 


UBA, ,UBA,, SESE 


J94uUao 393184 WOry YUeWOBTdSTp s,UBA —- 


LEO 


LAVYOULV 


pasa Sogn Kot Snopoaiaalh a'nev + way ‘iaw* av 59222” oe 


ot " Poaqeer We aoiziece & "stato71A < - + 2.6 @ = 
totus! 
o? sosyers \w aaibteds a Fists = ee Sr 
mer 


nearest of acl ht000 


i . 
[ne 7: = : ew 


C. BOMB DROP SIMULATION PHASE 

Prior to simulating the bomb drop, the pilot's release 
conditions, altitude, dive angle, and velocity, must be de- 
termined. , 
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(See Figure 6). 

‘At this point, the initial velocity and drop angle of 
the bomb are determined taking into account the ejection ve- 
locities (components perpendicular and parallel to the flight 
path) of the bomb rack canisters. 


Thus: 


VE3B 
+ ArcTan (JoLtVEDB? 


./(Vel + VE2B)2 + VE3B2" 


where: VE2B - the ejection velocity parallel to the air- 
craft's flight path. 


Drop Angle (a) 


Ul 


Velocity (Vel) 


VE3B - the ejection velocity perpendicular to the 
aircraft's flight path. 


In the case of a dual drag bomb (a projectile whose aero- 
dynamic characteristics change in flight), fin opening times 
must be considered. The purpose of a fin opening time is te 
determine when the bomb's fins open, thus altering the drag 


characteristics of the projectile. 
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Figure 6. Aircraft's Release Angle x and 
Bomb's Release Angle e 
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In the bomb drop simulation, it is assumed that the 
bomb falls straight down the aircraft's run-in-line. Since 
the bomb has no independent propulsion units and ballastic 
dispersion is not introduced in the theoretical model, this 
is a valid assumption. 

1. Differential Equations 

The mathematical model can be described in vector 

form as follows. Let x and y be the horizontal and vertical 


[f 
components respectively. From Newton's Equation cf Motion: 


a Bis Will 
eto coe val Ve Gp? 
ay We -gm - ||v,|| (Vv. * C,) 
dt t t its b 
where: 
Ve = total velocity 
m - mass of the projectile 
g¢ - force due to the earth's gravitation 


(V,*C, )- sum of the forces acting on the body, where 
dreg is proportional to the square of the 
velocity 

The drag coefficient (C,) is computed for each type 
of bomb and it depends on the speed at which the projectile 
is travelling, the cross-sectional area of the bomb, and the 


air mass density at a given altitude. 


In the component form, the above equation becomes: 
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where x and y represent horizontal and vertical displacements 
respectively, and where We and Ve are horizontal and vertical 
velocity components. 
a. Air Density Computation 
dp 


Given the hydrostatic equation AB (1) 


and the equation of state p oRT (2) 
‘where: p = pressure 

p = density of air 

Z = altitude 

g = gravitational force (32 ft/sec?) 

T = temperature 

R = gas constant 


Solving equation 2 for p and substituting into 


equation 1: 


dp Pa 

dz RT 

DS apre o Bs 

p dz RT 

l(a 0) ee 
Zz RT 


Take the integral thru some altitude Z: 


dz RT 


Z Z 
=0 z=0 


Assume T is constant. 


In p(z) - In p(0O) = - ae 
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Assuming g/RT to be constant, we reduce to 
—Z 
ie) Po 
where: H = Ba 
g 
Po density at z=0. 


Thus when z = H, 


o(H) = p(0) et 
Si ee GO) 
PCH Ss ae 
(Gy) ae 

(0) e 


Thus, we must choose H such that 


e(H) _ 1 
KO) 2 


e(0) = 0.00237691 slugs/ft* (Density of air at sea level.) 
Solving for e, 

e(H) = 0.0008744 slugs/ft? 
Solving for H from a fourth order density equation, 

H = 30500 ft. 


Therefore, 


305 
p(z) = Py © paeo 
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Using this formula, a divided difference density 
table was computed with the altitude as the independent vari- 
able. 

Similar divided difference tables were constructed 
for the speed of sound based on altitude and.the drag coef- 
ficient based on the Mach number of the aircraft. 

2. The Integration Scheme 
The system of differential equations was solved by a 
standard fourth order Runge-Kutta integration scheme which 
was selected for its computational efficiency. Let the four 


differential equations in the system be represented by: 
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The application of the Runge-Kutta scheme gives, for 
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The trajectory simulation utilizes this integration 
process to time step ene projectile towards the target. 
After every At incrementation, the altitude (y position) is 
compared with the target altitude. If the simulated bomb is 
above target altitude, the faveeemeaiiom a6 refreshed and up- 
dated. The present position and velocity components are 
saved for possible interpolation. If, because of a large 
time interval, the bomb travels Set eae altitude, a 
linear interpolation is used to reduce the time interval. 
The previously saved position and velocity components are 
restored and the integration is performed using the interpo- 
lated At. Impact occurs when the projectile altitude is 
within 1.2 inches of target altitude. At impact time, the 


values x, jy, vo and vy are retained for the next phase. 


D. IMPACT PHASE 

The impact point must be computed in two frames of refer- 
ence in order to be of benefit to the pilot. It must be cal- 
culated with respect to the target's run-in-line and the 
aircraft's run-in-line. 

Thus two run-in-line angles must be computed, an uncor- 
rected angle, ve a ArcTan ) for the aircraft's impact 
point, and a corrected yY due to the velocity component of the 


effective wind EW = (Xwind, Ywind) 


e re 2 
Y = ArcTan (wand for the target's impact points. 
Y + Ywind 


(See Figure 7). 
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Figure 7, Aircraft's Run-In-Line Angle % 
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Using DR as the down range travel of the bomb and Time 
as the total time for the bemb to reach the impact point 
(See Figure 8), the impact points are calculated as follows: 

1) the target's impact coordinates, 


us = @e Sy " a er 
1. DR Sin y Xa/e Xwind*Time 


I 


+ DR * Cos y - nue - Ywind*Time 


Je 
2) the aircraft's impact coordinates, 
= * = * i 
T(a/c) 1. Cos Vy Be Sin Vy 
= * ee * 
1 (a/e) I. Cos Yq he Cos Y4 


where: DR - down range travel of the bomb. 
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III. FORTRAN SIMULATION 


' A. BOMB DROP SIMULATION 

The initial algorithms for simulating the TRAJ routine 
were supplied by Lee Thomson of the Naval Weapons Center, 
China Lake, California. 

Euler's integration scheme was used in- updating velocity 
and position components. The time interval set set at .01 
seconds for accuracy and thus the number of calls made to 
the integration routine (EINT) was of eine order of 10%. Due 
to processing time restrictions, the Runge-Kutta method was 
selected as an integration scheme. The number of calls to 
EINT were hence reduced by a magnitude of one hundred, with 
results differing in the fifth decimal place. 

The TRAJ routine was tested using several release condi- 
tions to demonstrate the model flexibility in simulating 
pond loft (nuclear weapons), level flight release, and dive 
bombing conditions. Results compared favorably to published 
results with mean error rates of 1.15% for time of fall and 
.74% for down range travel. 

Choosing a proper time’ interval for the Runge-Kutta meth- 
od was conducted by varying the time interval from one to 
five seconds. The emphasis, at this point, was the number 
of calls to EINT for the projectile to reach target altitude. 
The results by altering At, showed changes in the sixth deci- 
mal place, and a substantial decrease in the number of calls 


to EINT. However, beyond a time interval of three seconds, 
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there was only a slight decrease in calls. This was due to 
the time interval interpolation necessary for the projectile 
to reach the target. The savings made in reduction of inte- 
gration calls were offset by increases in the interpolation 
section. The results of these test runs appear in Appendix 


a" 


B. OVERALL SIMULATION 

‘Divisions are the most time-consuming operations on the 
.-micro-computer. The number of divisions in the overall sys- 
tem was therefore greatly reduced by replacing constants 
used in the division operations by their reciprocals. This 
was done to replace the division process by the less time 
consuming multiplication process. 

The divided difference table look-up routine was imple- 
mented to save in the number of arithmetic operations. The 
routine is used to obtain values for the air mass density, 
the speed of sound, and the drag coefficients necessary in 
the TRAJ routine. The initial call to TABLE searches the 
table returning the proper value, and stroes the position in 
the table where the value was found. Subsequently, this 
position is used as the first estimate for table location on 
the next call for corresponding entries. Because the values 
of independent variables change by a small amount at-each 
iteration, the desired table entry is in the immediate vicin- 
ity of the last look-up. This method of computing functions 
trades execution time with memory space and allows the total 
calculation to be carried out in a matter of seconds of com- 


puter time. 
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IV. PL/M IMPLEMENTATION 


To implement the model on the micro-computer, a floating 
point arithmetic package was written to anit alin a four hexa- 
decimal digit accuracy during calculations. The floating 
point variable therefore occupies three, eight bit words; 
two words for the variable -mantissa.and one word for the man- 
tissa sign bit and the variable exponent. The mantissa of 
the variable is viewed be a normalized binary number with 
the left most digit always one. (See Figure 9). 

Tn order to manipulate these variables, an arithmetic 
subroutine library had to be built. At present, there exist 
routines that add and subtract, multiply, divide and take 

the square roots of these floating point variables. 
: The add routine compares the exponents of the inputs. 
After performing the proper mantissa shift, making the vari- 
able exponents equal, a full sixteen bit add operation is 
done with adjusted mantissa. To subtract, the mantissa sign 
bit is switched. | 

The multiply routine divides each of the input mantissas 
into four hexadecimal digits. The hex digits are then appro- 
priately multiplied using a table lock-up scheme for the hex 
product. The result is an eight hex digit number of which 
only the four most significant digits are used for the proda- 
uct mantissa. The fifth digit is used for rounding. The 
exponents are then added and the proper sign of the product 


is masked in the result. 
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The divide routine performs a full sixteen bit subtract 
on the input's mantissa. If the carry bit resulting from 
the subtraction is on, one is added to the quotient mantissa. 
The quotient mantissa is then left-shifted, one bit. This 
procedure is repeated sixteen times after which the original 
input exponents are subtracted and the proper sign is in- 
serted into the result. 

The square root routine utilizes Newton's Iteration for 
approximating square roots. At present, interrupts have not 
been implemented for negative inputs. | 

The actual PL/M version of the simulation program was 
translated directly from the working FORTRAN version. All 
associated flow charts, process graphs, program listings, 
and sample output can be found in the appendix. 

The PL/M routines have been compiled, however due to 
time constraints, none have been tested. Routines which 
still must be written for full system implementation are Sine, 
‘Cosine, and Arctangent trigonometric routines, and input/ 
output procedures. The I/O procedures are important in that 
they will act as the operator interface with the micro-com- 
puter during the actual simulation. 

All associated hardware interfaces are in building stages 
at the Naval Weapons Center in China Lake under the direction 


of Lt. Alfred Pease. 
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V. CONCLUSIONS 


The purpose of the thesis was to help develop a working 
system for practice bomb drop simulation and to demonstrate 
that micro-computers can carry out significant computational 
tasks. The reader will note that the hardware cost of the 
MCS 8 micro-computer to be used in the target system would 
be less than $1000. 

Because of the decreasing hardware costs of the micro- 
computer and significant increases in the processing speeds 
which the micro-computer industry is presently realizing, 
new systems which are equivalent to presently available 
avionics systems can be developed at a hardware cost which 
is nearly one hundred times smaller. The ee develop- 
ment costs are also likely to be smaller because the system 
can be subdivided into smaller, independently operating sub- 
systems, each of which perform a well defined task. 

As an example of future expansion, consider the No-Drop 
Bomb Simulation in conjunction with a release prediction sys- 
tem. Two micro-computers would be interfaced making it a 
multiprocessing system. One micro-computer would track the 
aircraft to the release position while the other micro-com- 
puter would perform a drop Simulation at the various track- 
ing positions. Thus, the pilot would constantly be informed 
of his position relative to the optimum release conditions. 


Systems of this nature are presently in operation, however 
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equipment costs are extremely high. The small cost of the 
micro-computer would make on board aircraft implementation 


possible. 
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APPENDIX A PROGRAM FLOW CHARTS AND PROCESS GRAPHS 


Symbols Used In Flow Charts 


Start/Stop 


Poe Input/Output 


Card Input 


Processing Symbol 
(Arithmetic Operation) 


Conditional Symbol 


Subroutine Call 
(ENTRY) - Entry Point 


Return Statement 


Program Entry Point (label) 
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Symbols Used in Process Graphs 
input output 


1, Logical Gates: 


AND Gate 
e 
jle> 


OR Gate (inclusive) 


XOR Gate (exclusive) 


+ Sue Shift Right I Bits 
a v.| shift Left I Bits 
ba 
See Shift Carry Bit into Word I Bit(s) Right 
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3. Bite Operations (8 bits): 


| ae | ADD 


x | MULTIPLY 
ieee OP. 
carry bit output 
| - | SUBTRACT 


ze DIVIDE 


>» ASSIGNMENT 


Cavinelont | eted 2G 
! 


low trartexe) 1400 0X 


— 


_ 


evi. I raut*ne tia 
oo I etal? esiae 
ie ie) a3 id voted: #7ide 


e\etle &) anuifereg 


QA 


[Pe Oe Ep 


4, Address Operations (16 bits): 


; Full 16 Bit ADD 
Ub 


Ee: | ee 


Full 16 Bit SUBTRACT 


o 


Returns 16 Bit Address of Variable 


5. Conversion Operations: 


Converts Byte to Address Word (upper 8 bits 4 0) 


HICH 


[2] [2 


Converts Leftmost 8 Bits of Address Word to 
Byte Word. 


[ noes 


Converts Rightmost 8 Bits of Address Word to 
Byte Word 


E 


6. Conditonal: 


[” false part 
—— CONDITION 


L, true part 


7. Branching: 


(5) Destination 
(4 Entry Point 
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8, Floating Point Operations: 


ADD ADD (Subtract) Routine 


MULT MULTIPLY Routine 


DIVIDE Routine 


oOo 
- 
lz | A] 


i) 
cee 
By 


SQUARE ROOT Routine 


9. Subroutine Calls 


TAS | Two Hexidecimal Table Look-Up (MULTIPLY Routine) 
Ines! Floating Point Table Look-Up Routine 

heaay Routine to Calculate First Derivate of Input 
pos 


Runge-Kutta Integration Routine 
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APPENDIX B PROGRAM DEFINITION OF VARIABLES 


SYMBOL 
PL/M if 
FORTRAN | varied 


A(T) RAO 


AA F1ISUM1 
F1ISUM2 
FISUM3 
AAA F2SUM1 
F2SUM2 
F2SUM3 


ACOEF 


ALPHA 


AMACH 


BAR(L) BARI 


BAR2 
BAR3 
CAT 
DCD 
DEG 
DELT 


Jo58 


RADIANS 


RADIANS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


RADIANS 


DIMENSIONLESS 


JME 
RADIANS 
RADIANS 
FT 


DIMENSIONLESS 


DEGS/RADS 


SEC 


DEFINITION 


Coefficients calculated 
using Legendre Polynomials 
for a Least Squares Posi- 
tion and Velocity Fit. 

I=1 - AO 
I=2 - Al 
I=3 - A2 for each r,@,9¢ 


Summation of least signifi- 
cant evaluated polynomials 
in Least Squares Fit. 


Summation of most signifi- 
cant evaluated polynomials 
in Least Squares Fit. 


Drag coefficient vaiues in 
inputed Mach vs. Drag 
table(s). 


Pilot's bomb release dive 
angle. 


Mach number values in in- 
puted Mach vs. Drag 
table(s). 


Base values for least recent 
entries in the position 
history registers. 


Target altitude 

Drag coefficient divided 
difference values in in- 
puted Mach vs. Drag 
table(s). 


Radians to degrees conver- 
sion constant 


Total fin opening time for 
dual drag bombs. 
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SYMBOL 
PL/M if 


FORTRAN |} varied DEFINITION 


DELT1 


DELT2 


DELTA 


DEN 


DIFF 


DR 


DRHO 


DT 


DTO 


DVSB 


DY(I) 


ENT 


ENTRY 


DY1 
DY2 
DY3 
DY4 


SEC 


SEC 


1/SEC 


VARIABLE 


FT 


SLUGS/FT® 


SEC 


SEC/FT 


FT/SEC? 
FT/SEC 
FT/SEC, 
FT/SEC 


VARIABLE 


VARIABLE 
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Fin opening time perpendic-— 
ular to flight path. 


Fin opening time parallel 
to flight path. 


Time interval for Runge- 
Kutta integration. 


Values used in calculating 
the first derivatives. 


Divided difference table 
values entered in Table 
look-up routine. 


Down range travel of the 
bomb at impact. 


Density divided difference 
values in inputed Alti- 
tude vs. Density table. 


Time interval for bomb simu- 
lation (subject to change 
during execution). 


Time interval for bomb simu- 
lation (fixed). 


Speed of sound (inverse) 
divided difference values 
in inputed Altitude vs. 
1/Speed of Sound table. 


First derivatives of bomb's 
Simulated position and 
velocity components. 


Independent table values 
entered in Table look-up 
routine. 


Independent variable for 
Table look-up routine. 
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SYMBOL 


’ FORTRAN 


GAMMA 


GAMMA1 


IPICKL 
IROCK 


JROCK 


Ki 


K2 


PL/M if 
varied 


RADIANS 


RADIANS 


Jae 


DIMENSIONLESS 


PICKL DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 
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DEFINITION 


Main: aircraft's run-in- 
line with respect to 
target's axis. 

Traj: release angle of the 
bomb adjusted for 
ejection velocities. 


Aircraft's wind corrected 
run-in-line with respect 
to target's axis. 


Altitude values inputed for 
Altitude vs. Density table 
Altitude vs. 1/Speed of 
Sound table. 


Pointer used in position 
history registers. 


Simulated bomb release sig- 
nal from the aircraft. 


=O single drag bomb 
=1 dual drag bomb 


Applicable to dual drag 
bombs. 
=1 first stage 
=2 second stage 


Pointer to release positions 
in the position history 
registers. 


Pointer for the Legendre 
Polynomials for Least 
Squares position and ve- 
locity fit. 


Pointer for the position 
history registers for 
Least Squares position 
and velocity fit. 
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BHA 


SYMBOL 


FORTRAN 


KGATE 


LKCD 


LKRHO 


LKVSB 


MACH 


MASS 


MAX 


MMACH 


NMACH 


PL/M if 
varied 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 


1/LB, 


DIMENSIONLESS 


DIMENSIONLESS 


DIMENSIONLESS 
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DEFINITION 


Applicable to dual drag 

bombs 

=1 ist stage (bomb @ 
release) 

=2 ist stage to fin 
opening 

=3 fin opening time 

=4 2nd stage 


Inputed pointer .to the 


table location of pre- 
viously found result in 
Table look-up routine. 


Pointer to table location 
of the previously found 
drag coefficient in the 
Table look-up routine. 


Pointer to table location 
of the previously found 
density in the Table 
look-up routine. 


Pointer to table location 
of the previously found 
speed of sound (inverse) 
value in Table look-up 
routine. 


Mach number of falling 
bomb (speed in relation 
to the velocity of sound 
at a given altitude). 


Inverse of the mass of the 
bomb. 


Inputed number of entries 
in associated table to be 
evaluated in Table look-up 
routine. 


of entries in Mach 
Drag table(s). 


Number 
vs. 


in the 
Drag 


Number of entries 
appropriate Mach vs. 
table. 
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SYMBOL 


PL/M if 

FORTRAN | varied 

NUM 

O(I,16) 01(16) 
02(16) 
03(16) 

POS(L) POS1 
POS2 
POS3 

Qs 

RAD 

RHOT 

Ss 

S1 V1ISUM1 
VISUM2 
V1iSUM3 

82 V2SUM1 
V2SUM2 
V2SUM3 


DIMENSIONLESS 


DIMENSIONLESS 


Jab 
RADIANS 
RADIANS 


SLUGS / SEC 


RADS/DEGS 


SLUGS /FT° 


FT? 


FT 
RADIANS 
RADIANS 


sae 
RADIANS 
RADIANS 
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DEFINITION 


Inputed pointer to appro- 
priate Mach vs. Drag 
table in Table look-up 
routine. 


Evaluated Legendre Poly- 
nomials for 


AO I=1 
Al T=2 
A2 T=3 


Polar release conditions 
‘resulting from Least 
Squares position fit. 


Temporary storage value used 
in calculating the first 
derivatives. 


Degrees to Radians conver- 
sion constant. 


Density values in inputed 
Altitude vs. Density 
table. 


Cross-sectional area of the 
bomb. 


Summation of the least sig- 
nificant evaluated poly- 
nomials multiplied by the 
corresponding values in 
the position history reg- 
ister in the Least Squares 
Fit. 


Summation of the most sig- 
nificant evaluated poly- 
nomials multiplied by the 
corresponding values in 
the position history reg- 
ister in the Least Squares 
Fit. 
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SYMBOL 
PL/M if 
FORTRAN | varied 


SO(1) SO1 


S02 
S03 


SY(1) SY1 
SY2 
SY3 
SY4 


TILDA(L) TILDA1 
TILDA2 
TILDA3 
TIME 


TRK(L,16) TRKi 
TRK2 
TRK3 
VAL 


VALUE 


VE2B 


VE3B 


DIMENSIONLESS 


FT/SEC 


1Ay 
RADIANS 
RADIANS 
SEC 


IAL 


RADIANS 


RADIANS 


VARIABLE 


VARIABLE 


FT/SEC 


FT/SEC 


81 


‘ 


DEFINITION 


Sum of the squares of the 
evaluated Legendre poly- 
nomials for 


AO T=1 
Al I=2 
A2 I=3 


Temporary storage of bomb's 
simulated positions and 
velecities. 
I=1 horizontal velocity 
I=2 vertical position 
I=3 horizontal position 
I=4 vertical velocity 


Inverse of the tracking 
time interval. 


Base values for most recent 
entries in the position 
history registers. 


Time for bomb to reach im- 
pact. 


Position history registers 
containing the 16 most 
recent r,9,¢ radar output 
values. Actual values 
contained are the differ- 
ences between the actual 
radar output and the most 
recent base value. 


Dependent table values en- 
tered in the Table look- 
up routine. 


Dependent variable (output) 
for the Table look-up 
routine. 


Ejection velocity parallel 
to the flight path. 


Ejection velocity perpendic- 
ular to the flight path. 


SYMBOL 


| PL/M if 

FORTRAN | varied DEFINITION 

VEL FT/SEC Total velocity of the bomb 
during the trajectory 
Simulation. 

VEL(I) VEL1 FT/SEC Polar release velccities 

VEL2 RADIANS/SEC resulting from the Least 
VEL3 RADIANS/SEC Squares velocity fit. 

VSBT SEC/FT Speed of Sound (inverse) 
values in inputed Alti- 
tude vs. 1/Speed of Sound 
table. 

VTY FT/SEC Total velocity of the air- 
craft at bomb release. 

XAC FT Aircraft's X position com- 
ponent at release, with 
respect to target center. 

XDOT FT/SEC Aircraft's X velocity com- 
ponent at release. 

XIX IMPX FT X impact position with re- 
spect to target's run-in- 
line. 

XIXAC IMPXAC FT X impact position with re- 
spect to aircraft's run- 
in-line. 

XIY IMPY FT Y impact position with re- 
spect to target's run-in- 
line. 

XIYAC IMPYAC FT Y impact position with re- 
spect to aircraft's run- 
in-line. 

XV FT X distance from van's axis 
to target center. 

XWIND FT/SEC X wind component. 

XX FT/SEC Temporary storage used in 
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calculating aircraft's 
release dive angle. 
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SYMBOL 
| PL/M if 


FORTRAN | varied DEFINITION 
YAC Aircraft's Y position com- 
ponent at release, from 
target center. 
YDOT FT/SEC Aircraft's Y velocity com- 
ponent at release. 
YM1(1) YM11 FT/SEC First estimate in Runge- 
YM12 Kutta Integration. 
YM13 
YM14 
YM2(1) YM21 . Second estimate in Runge- 
YM22 Kutta Integration. 
YM23 
YM24 
YM3(1) YM31 Third estimate in Runge- 
YM32 Kutta Integration. 
YM33 
YM34 
YM4(1) YM41 Fourth estimate in Runge- 
YM42 Kutta Integration. 
YM43 
YM44 
YT(1L) YT1 FT/SEC Updated estimate in Runge- 
YT2 FT Kutta integration, used 
YT3 FT to find associated first 
YT4 FT/SEC derivatives. 
YV FT Y distance from van's axis 
to target center. 
YY(1) yaal FT/SEC Bomb's simulated position 
YY2 FT and velocity components. 
YY3 Jak I=1 horizontal velocity 
YY4 FT/SEC I=2 vertical position 
I=3 horizontal position 
I=4 vertical velocity 
ZAC FT Aircraft's Z position com- 
ponent at release, from 
target center. 
ZDOT FT/SEC AircraLit"s Z velocaty com= 
ponent at release. 
ZV FT Z distance from van's axis 


to target center. 
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